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Abstract 

In this paper an application of Newman-Janis algorithm in spherical sym- 
metric metrics with the functions M{u, r) and e(n, r) has been discussed. After 
the transformation of the metric via this algorithm, these two functions M{u, r) 
and e{u,r) will be of the three variables u,r,9. With these functions of three 
variables, all the Newman-Penrose (NP) spin coefficients, the Ricci as well as 
the Weyl scalars have been calculated from the Cartan's structure equations. 
From these NP quantities, a class of rotating solutions of Einstein's field equa- 
tions can be obtained. These solutions include (a) Kerr-Newman, (b) rotating 
Vaidya solution, (c) rotating Vaidya-Bonnor solution, (d) rotating Husain's 
solution, (e) rotating Wang-Wu solutions. It is found that the technique de- 
veloped by Wang and Wu can be used to generate new embedded solutions, 
that the Kerr-Newman solution can be combined smoothly with the rotating 
Vaidya solution to generate Kerr-Newman- Vaidya solution, and similarly, Kerr- 
Newman- Vaidya-Bonnor solution of the field equations. It has also shown that 
the embedded universes like Kerr-Newman de Sitter, rotating Vaidya-Bonnor- 
de Sitter, Kerr-Newman- Vaidya-de Sitter can be derived from the general so- 
lutions with Wang-Wu function. All rotating embedded solutions derived here 
can be written in Kerr-Schild forms, showing the extension of Xanthopoulos's 
theorem. It is also found that all the rotating solutions admit non-perfect fluids. 

PACS number : 0420, 0420J, 0430, 0440N 
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1 Introduction 



In an earlier paper [1] it is shown that Hawking's radiation [2] could be expressed in 
classical spacetime metrics, by considering the charge e to be function of the radial 
coordinate r of Reissner-Nordstrom as well as Kerr- Newman black holes. Since these 
two black holes describe the 'stationary' metrics, it has intended to search for 'non- 
stationary' rotating metrics in order to incorporate relativistic aspect of Hawking's 
radiation in general relativity. The non-rotating Vaidya metric is a non-stationary 
generalization of Schwarzschild vacuum solution, describing the gravitational field of a 
null radiating star. Many attempts have been made to generate non-stationary rotat- 
ing metrics which describe the rotating external gravitational field of radiating bodies. 
To mention with, Vaidya and Patel [3] obtained a non-stationary rotating metric with 
mass M — —m{l + uip/a?) — (6^/4a^) cos^ ^}^/^ having minus sign, where m, a, h 
are constant and u is the retarded time coordinate. They claimed that their metric 
recovers the Kerr metric when 6 = 0. However, it is well known that the Kerr metric 
has mass without negative sign. Carmeli and Kaye [4] have shown, by considering the 
mass M of Kerr metric directly as function of coordinate u, that the Kerr metric can 
be made a non- stationary rotating metric, called variable mass Kerr metric. How- 
ever, Herrera and Martinez [5], Herrera et. al. [6] have discussed the interpretation of 
variable mass Kerr metric of Carmeli and Kaye. Gonzalez et.al.[7] have presented a 
non-stationary generalization of the Kerr-Newman metric, by allowing the three pa- 
rameters a, m and e to be functions of coordinate u and shown that the variable a{u) 
does not represent the rotating electromagnetic field of the Einstein- Maxwell equa- 
tions and concluded that to take the parameters of a metric as functions of u does 
not generalize the solutions enough. Mallett [8] applied Newman- Janis algorithm to 
the Reissner-Nordstrom-de Sitter 'seed' solution to derive a Kerr-Newman-de Sitter 
solution and afterward he considered the mass and the charged to be the functions 
of the retarded time coordinate u in order to get non-stationary charged radiating 
metric. Xu [9] has discussed the nature of the field equations of Mallett's solution. 
Jing and Wang [10] have considered the mass M{u) and the charge e{u) unchanged 
after the application of Newman- Janis algorithm to the non- rotating Vaidya-Bonnor 
'seed' solution with mass M{u) and charge e{u). In fact, after the application of the 
algorithm, the mass M{u) and charge e{u) should be functions M{u, 9), e{u, 9) of two 
variables u and 9. 

Here we employ the Newman- Janis algorithm [11] to generate rotating non-stationary 
metrics from the spherically symmetric 'seed' metric with the functions M(m, r) and 
e{u, r), where u and r are the coordinates of the spacetime geometry. Newman- Janis 
algorithm [11] is a complex coordinate transformation, which has been introduced by 
Newman and Janis to obtain Kerr metric, a rotating Schwarzschild vacuum solu- 
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tion of Einstein's field equations from the non-rotating Schwarzschild 'seed' solution. 
In another paper Newman et.al. [12] again applied the same transformation to the 
non-rotating charged Reissner- Nordstrom solution to get rotating charged Reissner- 
Nordstrom solution, which is now commonly known as Kerr-Newman black hole solu- 
tion in General Relativity. So this complex coordinate transformation can be used to 
derive rotating solutions from the non-rotating 'seed' solutions of Einstein's equations 
of spherical symmetric metrics. Herrera and Jimenez [13] applied the same Newman- 
Janis algorithm to an interior non-rotating spherically symmetric seed metric and the 
resulting rotating interior was tried to match with the exterior Kerr metric on the 
boundary of the source. Drake and TuroUa [14] have generated a class of metrics as 
possible sources for the Kerr metric by applying the same algorithm to any static 
spherically symmetric 'seed' metric. Drake and Szekeres [15] have shown the unique- 
ness of this algorithm in generating the Kerr-Newman metric and proved that the 
only electrovac Petrov type D spacetime generated by the algorithm with a vanishing 
Ricci scalar A is the Kerr-Newman space-time. Yazadjiev [16] has also shown that 
Sen's rotating dilation-axiom black-hole solution [17] can be derived from the static 
spherically symmetric dilation black hole solution via this algorithm too. 
The purposes of this paper are 

1. to apply the Newman- Janis algorithm to the spherical symmetric 'seed' metric 
with the functions M{u, r) and e{u, r) of two variables u, r, 

2. to calculate all the Newman-Penrose (NP) spin coefficients, the Ricci as well as 
the Weyl scalars [18] in general, 

3. to give examples of rotating solutions, published and unpublished, of Einstein's 
field equations from these NP quantities. 

The spherically symmetric metric with the functions of two variables u, r has been 
transformed via Newman- Janis algorithm [11] to get rotating metrics. After the 
transformation, M and e will be functions of three variables u, r, 9. Then we calculate 
all the Newman-Penrose (NP) spin coefficients, the Ricci as well as the Weyl scalars in 
general. Accordingly, the Einstein's tensors as well as the energy momentum tensors 
(EMT) of the matter fields have been presented in terms of complex null tetrad 
vectors. From this EMT one observes the description of having two fiuids system in 
the field equations. To visualize the two fluid system we rewrite the Einstein tensors 
in terms of one unit time-like and three unit space-like vectors constructed from the 
complex null vectors. 

Thus, we can generate rotating solutions mentioned in the abstract above from 
the Ricci as well as the Weyl scalars of the transformed metric. Consequently, some 
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of the results are cited for ready reference in the form of theorems based on rotating 
solutions discussed here. 

Theorem 1 // is the Kerr-Newman solution of Einstein's field equations and 

£a is geodesic, shear free, rotating and expanding null vector and one of the double 
repeated principal null directions of the Weyl tensor of g™ , then g™^ = g™ + 
2Q{u,r,6) iaib will he a rotating Kerr-Newman- Vaidy a solution with Q{u,r,6) — 
—rf{u)R~'^, where f{u) is the mass function of rotating Vaidya solution. 

Theorem 2 // g^^ is the rotating de Sitter solution of Einstein's field equations 
and ia is geodesic, shear free, rotating and expanding null vector and one of the 
double repeated principal null directions of the Weyl tensor of g'^f, then g™'^^ — 
g^^ + 2Q{r, 6) iah will be a Kerr-Newman-de Sitter solution with Q{r, 6) = —{rm — 
e^/2)i?~^, where m and e are constant and represent the mass and the charge of 
Kerr-Newman black hole. 

Theorem 3 All rotating stationary spherically symmetric solutions based on Newman- 
Janis algorithm are Petrov type D, whose one of the repeated null vectors, £a is 
geodesic, shear free, expanding as well as non-zero twist. 

Theorem 4 All rotating non-stationary spherically symmetric solutions, derivable 
from the application of Newman-Janis algorithm and possessing a geodesic, shear 
free, expanding and rotating null vector £a, are algebraically special in the Petrov 
classification. 

The Kerr-Schild ansatz of Theorem 1 can also be written in another form as 
5,KNV ^ _^ 2Q{r,9)£ah with Q{r,9) = -{rm - e^/2)R-'^, where m and e are 
constant and represent the mass and the charge of Kerr-Newman black hole and g^^^ 
is a rotating Vaidya solution of Einstein's field equations and is geodesic, shear 
free, rotating and expanding null vector of gf^. This ansatz can be interpreted as 
the Kerr-Newman black hole embedded into the rotating Vaidya null radiating back- 
ground, describing Kerr-Newman- Vaidya black hole. Similarly, theorem 2 states that 
Kerr-Newman black hole can also be embedded into the de Sitter cosmological back- 
ground, describing Kerr- Newman- dc Sitter black hole. Its alternative form will be 
the ansatz gf^'^^ = gf^ + 2Q{r,0)£aib with Q{r,0) = -(A*r76)i?-2, where A* is 
the cosmological constant. The extension of the theorem 2 in non-stationary version 
can be stated in the case of rotating Vaidya-Bonnor-de Sitter solution. Theorems 3 
and 4 follow from the stationary as well as non-stationary rotating solutions to be 
discussed in the next sections. 
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This paper is organized as follows: Section 2 presents a brief application of 
Newman-Janis algorithm to a spherically symmetric 'seed' metric with the functions 
M(u,r) and e(u,r). A general expressions of NP quantities with M{u,r,6) and 
e{u,r,6) are calculated from the Cartan's structure equations and cited for further 
use in section 3. The general properties of the rotating spherically symmetric metric 
is discussed after observing the nature of the NP quantities. Section 4 discusses the 
energy conditions of the energy momentum tensor for a time like observer with its 
four velocity vector. Section 5 gives examples of rotating stationary as well as non- 
stationary solutions from the NP quantities. In section 6 we introduce the Wang-Wu 
function in the rotating solutions, and it is shown that the general rotating solutions 
with Wang-Wu function can be used to generate rotating new embedded solutions 
like Kerr-Newman-de Sitter, Kerr-Newman-Vaidya, Kerr-Newman-Vaidya-de Sitter, 
rotating Vaidya-Bonnor-de Sitter. The conclusion of the paper is cited in section 7 
with suggestions and remarks of the solutions discussed in the earlier sections. 

The presentation of this paper is essentially based on the Newman-Penrose (NP) 
spin-coefficient formalism [18]. The NP quantities are calculated from Cartan's struc- 
ture equations written in NP formalism by Mcintosh and Hickman [19] in (-|-,-,-,-) 
signature. 



2 Newman-Janis algorithm 



For application of Newman-Janis algorithm, we start with a spherical symmetric 
'seed' metric written in the form 

ds^ = e^^du^ + 2dudr - r\de^ + siv?e d(j)'^), (2.1) 

where e'^'^ — l — 2M{u, r)/r + e'^{u, r)/r^ and the coordinate chosen are {x^, x^, x^, x^} 
—{u,r,9,(f)}. The u-coordinate is related to the retarded time in flat space-time. 
So u-constant surfaces are null cones open to the future. The r-constant is null 
coordinate. The 9 and are usual angle coordinates. The retarded time coordinate 
arc used to evaluate the radiating (or outgoing) energy momentum tensor around the 
astronomical body [9]. Here M and e are the functions of the retarded time coordinate 
u and the radial coordinate r. Initially, when M, e are constant, this metric provides 
the Reissner- Nordstrom solution and also when both M, e are functions of u, it 
becomes the non-rotating Vaidya-Bonnor solution. The contravariant components of 
the metric (2.1) are 
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a,h= 1,2,3,4. 
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These metric components can be expressed in terms of complex null tetrad [18] 



1 2M 

n'-Si--(l- — + ^)Sl (2.3) 

Then the metric tensor g"'^ of the line element (2.1) is expressed in these null tetrad 
vectors as 

g"-^ ^t-n!' + rft -rrfrf^ -nfm!'. (2.4) 

Here the null vectors and are real, and rrf" and m" are complex conjugates of 
each other. According to Newman and Janis [11], one can complexify the coordinate 
r and u by the following transformation 

r = r' — ia cos 9, u = u' + ia cos 9, 9 = 9', 4> = 4>'. (2.5) 

This complexification can only be done by considering r' and u' real. Thus, from 
(2.5) we have the relations 

dr = dr' + iasm9 d9, du = du' — iasm9 d9, d9 = d9', and dcj) = dcj) . (2.6) 

These relations can be utilized in (2.1) to obtain the rotating metrics: 

ds'^ = e^^du^ + 2dudr + 2a sin^ 9 {l-e^'^)dud<j) 

-2a sm'^9 dr dcf) - B?d9'^ - {R^ - a^ sm'^9 (e^"^ - 2)} sin^^ (2.7) 

where e^"^ = l-2rM{u,r,9)/R'^ + e^{u,r,9)/R^ and R^ = r'^ + a^cos^9. All the primes 
are being dropped for convenience of notation. According to Newman and Janis [11], 
we have also used a suitable substitution d9 = —i sm9 dcf). Then the covariant complex 
null tetrad vectors take the forms 

la^5l-as\u^9 5l 

1 1 



n, 



a 



- H{u, r, 9) 51 + 51-- H{u, r, 9) a sm'9 5^, 
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la sin e5l + R^5l + i{r^ + a^) sin ^ 5^}, 



rria 



V^R 



1= {ia sin eSl + R' 61 -i{r' + a') sin 96^}. 



where R = r + ia cos6. The null tetrad vectors chosen here are different from those 
chosen in [12], but are similar to those given in Chandrasekhar [20]. Now, after the 
transformation (2.5), the functions M and e must be of the three variables u,r,9, 
however the old ones had explicitly u and r dependence. That is. 



2rM{u,r,e) e'^{u,r,e) a^sin^O 



i?2 i?2 ' i?2 ■ 

Then the line element has the covariant components of the metric tensor Qah- 



(2.9) 



9ah 



( e^'^ 10 a^\^Q{\-e^-^) \ 

1 -asin^^ 



V a sin^^ (1 - e^"^) -a sin^^ -{R^ - o? wi^d{e^'^ - 2)}sin2^ j 



(2.10) 



a,b — 1,2,3,4. This completes the application of Newman- Janis algorithm to the 
spherically symmetric 'seed' metric (2.1). The usefulness of this transformed metric 
(2.7) will be discussed in the following sections. 



3 NP quantities for the rotating metric 

In this section we shall derive the general NP spin coefficients, the Ricci scalars and 
the Weyl scalars for the metric (2.7) and present the general properties of the metric 
after observing the conditions of these NP quantities. First, the basis one-form of the 
tetrad vectors (2.8) are given below: 

6*^ = 4 d'^"' ^ du- a sin^O dcf), 
= -madx''^^^{iasinedu + R^de-i{r^ + a'^)d(l)}, (3.1) 

9^ = -madx" ^ -ia sine du + R'^d0 + i{r'^ + a'^)d(f)}. 

y2 R 
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The intrinsic derivative operators for the metric (2.7) take the following forms: 



D = tda = dr, 



5 = m''da 



+ 



^ {iasinedu + de + ^^dA, 



(3.2) 



5 = rrf da = -^r= \ - ia sin ed^ + do — \. 



V2i? 



By taking the exterior derivative of basis one-forms (3.1), we calculate the spin co- 
efficients from the Cartan's equations of structure written in Newman-Penrose spin 
coefficients [19]: 



H{u, r, 



a — 



{2ai — Rcos9) 
2^2RRsm9 ' 

i a sin ^ 

TT = , , T 

y/2RR' 



/3 



cot^ 

i a sin 9 
^2R^ ' 



(3.3) 



7 = 



^2RR^ 



{r - M - r M^r + e e,^) -R - A' 



^2RR^ 



i a sin ^(r - e e,„) - (r - e e,^) 



where A* = — 2rM{u, r, 9) + + r, 9) and the function H{u, r, 9) is given in 
(2.9). From these NP spin coefficients we conclude that the transformed metric (2.7) 
with the functions M and e of three variables u, r, 9 possesses, in general, a geodesic 
{k — e — 0), shear free {a — 0), expanding {9 ^ 0) and rotating (a;*^ ^ 0) null vector 
la [20] where 

^^^-^(P + P) = ;^, (3.4) 
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1 2 a^cos^^ 



.*2 



Further we calculate the Weyl scalars: 



(3.5) 



^2 



1 



RRE? 



-RM + e^)+R {rM^r - ee,^) + i i? i? ( - 2M^ - rM^rr + ej, + ee^rr) 



V'3 = 



[4{i a sin 9 (rM^^ - ee,„) - (rMg - ee,^) } 



2^2RRR'^ 

+R {i a sin e{r - e e,„),^ - (r - e e,^),^}] , 



(3.6) 



1p4 



1 



2RRR'^ 
r a^sin^^ 

~RRR^ R^ 



a^sni^O {rM^u - ee^u),u + 2iasin 9 {rM^u - ee^u),e - {rM^e - ee,^),^ 



(rM„ 



ee . 



2rasin2^ - i^^cos^ 
2RRR'^ R^sm9 



(rM 



The non- vanishing of the Weyl scalars {1^2 7^ V's 7^ "^4 7^ 0) means that the metric 
(2.7) is an algebraically special in the Petrov classification. In the expression of 1^2 it 
is found in general that there is no differential terms of M{u, r, 6) and e{u, r, 6) with 
respect to u and 9. This leads that for a static rotating metric with the mass M(r) and 
the charge e(r), the spacetime metric will be a Petrov type D {ip2 7^ 0, -03 = -04 = 0), 
whose one of the repeated principal null vectors is a geodesic, shear free, expanding 
(3.4) and rotating (3.5) vector £a. 

The Ricci scalars of the metric (2.7) are obtained as follows: 

000 = 001 = 010 = 020 = 002 = 0, 



1 



^11 



4i?2 R2 



2 + 4r(r - ee ^) + - 2M^ - r M,^^ + ej. + ee,^^) 



012 = 



i a sin 6l|(i?Mu - 2ee,„) - (rM^ - ee,r),«-R} 



2^2 i?2 R2 

+{{RM^e-'2ee,e)-{rM,r-eer),eR} 







'21 



;asin^{(i?M„ -2ee,„) - (r - e e,^),„ i?} 



2^2 i?2 i?2 

+ {(:RM, - 2 e e,,) - (r M,, - e e,,),, r] 



(3.7) 
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022 = 



2r(r M„ - e e,„) - cot ^(rM,^ - e e,^) 



sin^6'(r M„ - e e,„),„ - (rMg - e e,^),^ 



A = (2M, + r M,, - e^^ - ee,,,) . 



Here we observe that the general expressions of and A do not involve any 
differential terms with respect to u and 9, although the functions M and e of the 
metric (2.7) are functions of three variables u,r,9. The vanishing of 0oo suggests 
the possibility that the transformed metric (2.7) does not possess a perfect fluid 
Tab = (p* +p)uaUb-pgab as 000 = 20ii = 022 = -K{p* + p) / 4:, A = K{3p - p*)/24 
with a time-hke vector = {£"■ + n")/y2 [21,22]. It is also worth mentioning that 
for a static rotating metric with M(r) and e(r), the Ricci scalars 0i2 and 022 will be 
vanished. 

Then the Einstein's tensor is computed from these Ricci scalars (3.7) as follows 

Gab = -2 022 4 ^6 - 4 011 + m(„m6)} -GAgab 

+4 012 i{a rrib) + 4 02i i(a rrib) , (3.8) 

where 2i(^a^b)=^a^b + ^a4- Foi' non-rotating fields (a = 0), the Ricci scalars 0i2, 
021 will vanish and this Einstein's tensor will reduce to that presented by Glass and 
Krisch [23]. Prom the Einstein's equations: 

Gab = Rab — -Rgab — —KT^b (3.9) 

wc obtain the null density /x*, the matter density p*, the pressure p as well as the 
rotation function u) as 

ir//* = 2022, = -2 012, Xp* = 20ii + 6A, i^p = 20ii-6A (3.10) 

where the Ricci scalars 0ii, 0i2, 022, A are given in (3.7). 

To have the two- rotating fluid description we shall introduce a time- like unit vector 
u"- and three unit space-like vectors v"-, w"-, z"- such that 

Ua = -^{(-a + ria), 
Va = ^(4 - ria), 
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(3.12) 



1 i 

Wa= -^{ma + ma), Za = --j^{ma-ma) (3.11) 

with the normahzation conditions UaU'^ = 1, VaV"' — WaW^ — z^z"" — —1. Then the 
exphcit forms of these unit vectors are as follows 

va = ^^{{\-\ii)?>\-bl-{\-\ //) asin^^5„^}, 
Wa = --^ I - a^sin ^ COS ^ (5^ + ri?^ 5^ + a (r^ + a^)sin ^cos 6^ 

Za = jar sin 6* 5^ + a cos 6 R^Sl — r (r^ + a^)sin^5^|, 
The metric tensor g^b can be expressed in these unit vectors 

gab = UaUb- Va Vf, -WaWb- Za Zb- 

Thus the Einstein's equations are written in two-fluid system 

Gab = -K[^* iah + P* (UaUb-VaVh) +p{WaWb + ZaZb) 

+ {U + UJ) {u(aWb) +V(aWb)} - i{uJ - uj) {u(a Zb) + V(a Zb)}], 

where //*, p*, p and u; are related with the Ricci scalars given in (3.7) as: 

^A** = -j^{MrM,u-eeu)-cote{rM,e-eee) 
sin^6l(r M„ - e e,„),„ - (rM,^ - e e,e),e} 



(3.13) 



(3.14) 



^P= Wmi""^ + 2r{rM,, - ee,,) - R\2M^r + rM,,, - - ee,,,)}, (3.15) 



i a sin 6'{(i?M„ - 2ee„) - (rM^ - ee,r),„-R} 



i?2 i?2 

+{ (i? Me - 2 e e,e) - (r M,, - e e,,) ,e 



The expression of null radiation density fi* involves the derivative of the functions 
M{u, r, 6) and e{u, r, 9) with respect to u and 9. Those of p* and p are with respect to 



11 



r only. However, the expression of the rotation function uj is involved the derivative of 
mass and the charge with respect to three variables r, 9. Prom the above equations, 
it is observed that the Einstein's tensor (3.14) of the rotating fluid reduce to those of 
Glass and Krisch [23] of non-rotating string fluid when a = e = and M = M{u,r). 
In the non-rotating Vaidya-typc radiation null fluid, the null density /j* takes the form 
fi* = —2Mu/Kr^. This shows that /i* is always negative, since dM/du is positive 
[24]. 

4 Stress-energy tensor and energy conditions 

Prom the Einstein tensor (3.8) and the relations (3.10) of Ricci scalars with //*, p*, p, 
u we can introduce the total energy momentum tensor (EMT) for a rotating fluid as 
follows: 

rp _ rp(n) rp{va) 

J- ah — J-ah ^ -'-ab 

= II* iaib + '^P* ^{arib) + 2pm(^Jnb) 

+2u}£(arnb) + 2cj£(^amb) (4.1) 

where the EMTs for the rotating null fluid as well as that of the rotating matter are 
given respectively below: 

r^6^ = //* 4 4 + u; £(a ffib) + u}£{a rub) (4.2) 

T^^^ = 2 (p* + p) £(a nfc) -pQab + UJ i(^a rUb) + CJ £(a "Ife) , (4.3) 

where uo is the complex conjugate of uo. When = initially, these EMTs is similar 
to those introduced by Husain [25] in the case of non-rotating fluid. 

In General Relativity the stress-energy tensor represents the matter that describes 
the gravitation in the space-time geometry through Einstein's field equations. Prom 
the conditions of Ricci scalar ^qo = 0, obtained above for the space-time geometry 
(2.7), it may be concluded that the stress-energy tensor given above does not, in 
general describe a perfect fluid, as for a perfect fluid the Ricci scalar </)oo = —K{p* + 
j9)/4 must not vanish. Hence, it will be interesting to study the nature of the energy 
conditions for rotating non-perfect fluid given in (4.1). When the rotation factor 
UJ vanishes in (4.2), this fluid may be thought of null radiation fluid of non-rotating 
Vaidya space-time. So we refer to this rotating {uj ^ 0) null radiation fluid as rotating 
Vaidya type radiating fluid, shortly rotating Vaidya fluid (4.2). 

As the Tab does not include the perfect fluid, it seems that the stress-energy tensor 
can represent the interaction of rotating Vaidya fluid with rotating non-perfect fluid. 
Since there is a coupling term of the rotation scalar a with DM/ du in the expression of 
u) appearing in Tab, the energy condition of this Tab satisfying Einstein's held equations 
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will be a new area to discuss in the classical General Relativity. For this purpose, 
we write the matter part of (4.1) in terms of time-like u"' as well as space-hke 
vectors v"', w"-, z"- given in (3.12) as 

Tab = 1^*^^+ {P* +P)iUaUb-VaVb) -pgab 

+ {llI + U) {U(^aWb) + V(aWb)} - i{uj - Uj) {U(^a Zb) + V(aZb)}, (4.4) 

and its trace is T = TabQ""^ = 2{p* — p), which is different from that of a perfect fluid. 
This trace will be vanished when p* = p. This means that the matter part T^™"* of the 
stress-energy tensor may be that of electromagnetic field whose trace is zero. [The 
stress-energy tensor for a non-rotating perfect fluid is T^^^^ = (p* +p)ua Ub —pQab with 
unit time- like vector Ua and trace T^^^^ — p* — 3p, which will be zero when p* = 3p] . 

To study the energy conditions of the energy-momentum tensor, we consider a 
time-like observer with its four-velocity vector [/„ [26,27] 

Ua = aUa + (3Va + ^Wa + fca, (4.5) 

where a, /5, 7 and 5 are arbitrary constants. The four-velocity vector [/„ is subjected 
to the condition that 

U^Ua = - - 72 - ^2 > 0. (4.6) 

Now, TabU"'U^ will represent the energy density as measured by the time-like observer 
with the unit tengent vector f/". Then the energy conditions [36] are the following: 

(a) Weak energy condition: The energy momentum tensor obeys the inequality 
TabW^U^ > for any timelike vector i,e., TabW'U^ > imphes that 

(i)y + P*>0, (ii) ±f^ + p*+p>0, (m)^+p*+p>0, 

(iv) {^ + p*+p)±{u; + cJ)>0, (v) {^ + p*+p)±i{u;-u;)>0. (4.7) 

(b) Strong energy condition: The Ricci tensor for Tab (4.4) satisfles the inequality 
RabWU^ > for aU time-hke vector C/«, i.e. TabUW^ > \T, which imphes that 

(i)y+P>0, (ii) ±^ + p*+p>0, (m)^+p*+p>0, 

{iv) (^ + p*+p)±(u; + cz;)>0, {v) (^ + p* + p) ± i (a; - cZJ) > 0. (4.8) 

(c) Dominant energy condition: For all future directed, time- like vector U"", TabU^ 
should be a future directed non-space like vector. This energy condition is equivalent 
to 

ii) pV* + p*2 - [f + /) > 0, 
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{ii) - /.V* + - + if + > 0, (iii) + -p'- 3{f + > 0, 



where 2f — oj + uJ and 2i g = 00 — oJ. These are the energy conditions satisfied by 
the energy-momentum tensor (4.4) in general. We find that the rotation function u 
is involved in all the three energy conditions. 

5 Rotating solutions recovered from the general 
solutions 

In the above section we present the full expressions of NP spin coefficients (3.3) the 
Weyl scalars (3.6) and the Ricci scalars (3.7) with arbitrary functions M and e of 
three coordinate variables u, r, 9. These NP quantities are so transparent that these 
quantities can explain the nature of solutions of Einstein's field equations discussed in 
this paper. For example, the NP spin coefficients (3.3) easily explain that the metric 
(2.7) admits a null vector which is geodesic, shear free, expanding (3.4) as well as 
rotating (3.5), In this section, with the help of NP quantities given in (3.3), (3.6), 
(3.7), we shall give known examples of rotating solutions like Kerr-Newman, rotating 
Vaidya [4] and rotating Vaidya-Bonnor [10]. In the next section we shall combine 
these solutions with other rotating solutions in order to derive new embedded rotating 
solutions. 

5.1 Kerr-Newman solution: e = M = constant, a 

When e = M = constant, a 7^ 0, the equation (3.10) reduces to the Kerr-Newman 
solution 



(tv) {fx*p* + p*' -p'- if + g')} ± 2f{p* - p) > 0, 



(4.9) 



(v) {p*p* + p*' - / - if + /)} ± 2gip* -p)>0, 



p* ^ou^O 



p =p = 



(5.1) 



and the only existing Weyl scalar is 



Then the total energy momentum tensor takes the form 



Tab = P* (4 rih + ria 4) + P {rnaTrib + rua rub) 
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= {e^/KR'^R'^){{ianb + naib) + imamb + mami,)}, (5.2) 



which is the EMT for non-null electromagnetic field with Maxwell scalar 

for Kerr- Newman solution. Here is the birth place of Kerr- Newman solution, origi- 
nally applied the Newman- Janis algorithm by Newman et. al. [12] to generate this 
well known rotating solution from the non-rotating Reissner-Nordstrom 'seed' solu- 
tion. The line element is 

ds^ = {l-(2rM -e^)R~^}du^ + 2dudr 

+2aR~'^{2rM - e^)sin^6' du dcf) - 2a sin^O dr d(p 

-R^de^ - {(r^ + a'f - A*a2 sin^^ ji^'^sin^^ #^ (5.4) 

where A* = r'^ — 2rM -\- a? -\- . The charged Kerr-Newman black hole has an external 
event horizon a,t — M + -\/(M^ — — e^) and an internal Cauchy horizon at 
r_ = M — a/(M^ — — e^). The stationary limit surface Quu > of the rotating 
black hole i. e. r = re{6) = M + v^(M^ — a^cos^^^ — e^) does not coincide with the 
event horizon at r+ thereby producing the crgosphcre. This stationary limit coincides 
with the event horizon at the poles 9 = and 6 = n [20]. Naturally, this solution 
includes Kerr (e = 0), Reissner-Nordstrom (a = 0, e 7^ 0) as well as Schwarzschild 
(a = e = 0) solutions. 

5.2 Rotating Vaidya solution: M = M{u), a 0. e = 
In this case the energy momentum tensor (4.2) takes 

Tab = 11* iah + i^^iarnb)+uJ£(amf,) (5.5) 
where the null density /x* and the rotation function cu in (3.15) become 

Kfi* = --^{2r2M,„ + aVsin2eM„„}, 



and the Weyl scalars are 



Kuj = --^^^^iasinOM^^, (5.6) 



, _M_ 
RRR' 



?^a_sing UrM^ + RMA, (5.7) 

2y/2RRR^^ ' 'J' ^ ^ 
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From the above we observe that uj, t('3. ip^ will vanish when a = 0, and the EMT 
will be that of the original non-rotational radiating Vaidya metric [28] with //* = 
—2M^u/Kr^. The line element of this rotating metric is 

ds'^ = {l-2rM{u)R'-^}du'^ + 2dudr 

+AarM{u)sm'^9 du dcf) - 2a sm^9 dr dcj) 

-R'de^ - {(r' + aY - A* a' sin^^ji^-^sin^^ dcj)'', (5.9) 

where A* = — 2rM{u) + o?. This metric represents a non- stationary rotating 
solution of Einstein's equations possessing an energy- momentum tensor (5.5) for a 
rotating null radiating fluid, and can describe a non-stationary rotating black hole 
if M{u) > a. The involvement of uu in the energy-momentum tensor (5.5) indicates 
that the null fluid is a rotating Vaidya null fluid. When M (u) = constant initially, 
this metric would reduce to rotating vacuum Kerr solution with vanishing /i* and cu 
in (5.6). 

Carmeli and Kaye [4] studied the metric (5.9) after considering the mass M of 
the Kerr solution as a function of coordinate u. That is why, they referred to the 
metric (5.9) as the variable-mass Kerr solution ( see also in [29,30]) and discussed 
the properties of the metric using the NP quantities. Carmeh [29] referred to these u; 
terms as residues of the black hole. However, we refer to the metric (5.9) as rotating 
Vaidya solution. In the next section this metric will be combining smoothly with the 
usual Kerr-Newman solution as rotating Kerr-Newman- Vaidya black hole. 

5.3 Rotating Vaidya-Bonnor solution: M = M{u), a ^ 0. e = e{u) 
In this case the energy momentum tensor takes 

Tab = Ai*44 + 2p* {£(anb) + mi^Jfib)] + 2a;£(amb) + 2a;£(amb) (5.10) 

where 

A** = ~ D2 P2 ^ ~ ^ ^.") + a^siv?d{r M^u - ee,u),u\-, 



UJ — 



-i a sin 9 r ^ , ^ ^ 1 
|i?M„-2ee,„|, 



^2KR'^ 
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and the Weyl scalars are 



V'3 = 5^^^^{4(rM,.-ee„)+7iM„}, (5.12) 

The hne element will be in the form 

+2ai?-2{2rAf(M) - e^(M)}sin2^ rfw rf0 - 2a sin^O dr del) 

-R^de^ - {(r' + 0^)2 - A*a2 sin^^ ji^'^sin^^ d^^, (5.13) 

where A* — r"^ — 2rM{u) + + e^{u). This solution will describe a black hole when 
M(u) > + e'^{u) and has r± = M(u)* ± ^{M2(ii) - - e'^(u)} as the roots of 
the equation A* = 0. So the rotating Vaidya-Bonnor solution has an external event 
horizon at r+ = M{u) + ^{A^Piii) — — e2(-u)} and an internal Cauchy horizon 
at r_ = M{u) — ^{M'^{u) — a^ — e'^{u)}. The non-stationary limit surface Quu > 
of the rotating black hole i.e. r = re{u,6) = M{u) + ^J{M'^{u) — a^cos^d — e^{u)} 
does not coincide with the event horizon at r+, thereby producing the ergosphere. 
The rotating Vaidya-Bonnor metric (5.13) can be written in Kerr-Schild form on the 
rotating Vaidya null radiating background as 

gl''^gl + 2Q{u,r,9)ej, (5.14) 

where 

Q(n,r,e) = '^. (5.15) 

Here, is the rotating Vaidya metric (5.9) and is geodesic, shear free, expanding 
and rotating null vector for both as well as g^,^ and given in (2.8). The Kerr- 
Schild form (5.14) may be interpreted as the existence of the electromagnetic field 
on the rotating Vaidya null radiating background. If we set M{u) and e{u) are both 
constant, this Kerr-Schild form may be that of Kerr-Newman black hole. That is, the 
Kerr-Newman solution itself has the Kerr-Schild form on the Kerr background with 
the same null vector £a (2.8). 

From this rotating Vaidya-Bonnor metric, we can clearly recover the following 
solutions: (i) rotating Vaidya metric (5.9) when e(-u) = 0, (ii) rotating charged Vaidya 
solution when e{u) becomes constant, (iii) the Kerr-Newman solution (5.4) when 
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M{u) — e{u) — constant and (iv) well-known non-rotating Vaidya-Bonnor metric 
[31] when a — 0. It is also noted that when e — a — 0, the null density of Vaidya 
radiating fluid takes the form fj,* = —2M^u/Kr'^. The non-rotating Vaidya null 
radiating metric is of type D in the Petrov classification of spacetime whose one of 
the repeated principal null vectors, is a geodesic, shear free, non-rotating with 
non-zero expansion [29], while the rotating one is of algebraically special with a null 
vector ^a which is geodesic, shear free, rotating as well as expanding. It is also noted 
that when e = a = 0, the energy-momentum tensor becomes that of the original 
non-rotational null-radiating Vaidya fluid with ji* — —IM^y^jKr^. It is noted that 
the metric (5.13) can be seen in [10]. Prom (5.4) and (5.13) it is observed that the 
rotating Vaidya-Bonnor solution is the non- stationary version of Kerr-Newman black 
holes. That is, the parameters M and e of Kerr-Newman solution are functions of 
retarded time coordinate u in rotating Vaidya-Bonnor metric. 

6 Rotating solutions with M = M(u, r), e{u, r, 6^) = 

In this section we shall discuss the rotating solutions by considering the function 
M{u, r) of two variables u, r only and e{u, r, 9) — 0, and derive new embedded rotating 
solutions. We shall express all the rotating embedded solutions in Kerr-Schild types 
of metrics in order to show them as solutions of Einstein's fleld equations. In this 
case the energy momentum tensor takes the form 

Tab = + 2p*^(ari6) 2pm(amfe) 2a;£(amfe) + 2a;£(amfe) (6.1) 

or in terms of unit vectors 



A** 4 (-b + P* {Ua Ub - Va Vb) + P {Wa Wb + Za Zb) 

+ {UJ + Uj) {U(a Wb) + V(a Wb)} - i{uj - Uj) {u(^a Zb) + V^a Zb)}, 



(6.2) 



where 




M 



KR^ i?2 




(6.3) 



The line element will take the form 
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+4:arM{u, r)R~'^ sva^O du d(f) - 2a sin^O dr d(f) 

-RHe'' - {(r^ + a^f - A*a'' sm'^O} R-\in^ dcf)'', (6.4) 

where R^ = r"^ + a^cos^O, A* = — 2rM{u, r) + and the Weyl scalars are 

1 ~R T 

^2 = ^ { -RM+—MJAr + 2iacose)--^'RMrr\, 

^ RRR^^ 6 ^6 'J' 

One can consider this rotating metric (6.4) along with the stress-energy momentum 
tensor (6.1) or (6.2) and the Weyl scalars as the extension of the non-rotating solutions 
discussed by Glass and Krisch [23] and Husain [25]. 

6.1 Rotating Husain's solution: M = M{u,r), a 

Husain [25] has imposed one condition in the equation of state of non-rotating null 
fluid that p — kp*^ and obtain the solution of the equation of state with k > 1/2 
b — 1. However, due the present of the rotating factor a in equation (6.3), one cannot 
be able to get the solution. So we put k — 1 and b — 1. Then, the equation to be 
solved takes a simple form 



Mr Mrr , x 

-y - (6-6) 

which gives the function M{u, r) 

M{u,r)^f{u)--g{u). (6.7) 
r 

It can be treated as rotating Husain's solution for p = kp* with k = 1. This ro- 
tating Husain's solution may be degenerated to the rotating Vaidya-Bonnor solution 
presented above if one puts g{u) = e^{u)/2 in (6.7). 

6.2 Rotating Wang-Wu solutions 

Wang and Wu [32] have expanded the function M{u,r) of (6.3) of the non-rotating 
space in the power of r 

M(«,r)= Qn{u)r'', (6.8) 

n=—oo 

where qn{u) are arbitrary functions of u. They consider the above sum as an integral 
when the 'spectrum' index n is continuous. In fact Wang and Wu technique is based 
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on a linear superposition that a linear superposition of mass function of particular 
solutions is also a solution of Einstein's field equations of non-rotating spacetime. 
Using the expression (6.8) in equations (6.3) we can generate rotating solutions with 
Wang-Wu functions as 



^ {2g„(«),„r"+i + aW^g„(«),„„r"}, 



2 J.2 +00 
±\ ±L ±L n=-oo 



p 



E ng„(M)r" ^1 + (n-l)|, 



KR^ ^fr'oo ' ^ 



i a sin 6 / — \ . s „ , ^ 

/2KR^R^ E {R-nR)qn{u),^T\ (6.9) 

V n=—oo 



Here one can observe that these rotating solutions with functions qn{u) include many 
known as well as unknown rotating solutions of Einstein's field equations in spherical 
symmetry as shown by Wang and Wu in non-rotating cases [32]. The functions 
qn{u) in (6.8) play a great role in generating new solutions whether rotating or non- 
rotating. Therefore, we will hereafter refer to these as Wang-Wu functions. Thus, 
rotating solutions can be derived from these solutions as follows. 

6.2.1 Rotating monopole solution 

If one chooses the functions qn{u) such that 

, , / (6/2), when n = l 

= I 0, when n ^ 1 ^^'^^^ 

where b is constant, then from (6.8) one can obtain 

M{u, r) = br/2, n* = cu = 0, 



^ r'^b ba^cos^O 
^ " KR^R^' KR^R^' 

with the energy momentum tensor 

Tab = 2 p* ni,)+2p m^Jfih) ■ 
The Weyl scalar takes the form 

br 

~ ~2RRR' 
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The rotating monopole line element will be of the following form 

ds^ = {l -br^R-^}du^ + 2dudr 

+2a b r^R~^ sin^O du dcj) — 2a sin^O dr dcp 

-R^de^ - {(r' + - A*a^ sin^O} R-^m^O d(t)\ (6.11) 

where R^ = r"^ + a^cos^O, and A* = — + a^. 

Because of the non- vanishing Weyl scalar il)2-, the rotating monopole solution is 
stationary Petrov type D whose one of the repeated principal directions is geodesic, 
shear free, expanding and rotating null vector The rotating monopole solution has 
the non-zero pressure p, which leads the difference between the rotating and the non- 
rotating monopole solutions. That is, when a = 0, one can obtain the non-rotating 
metric [32] with p = 0. To study this rotating monopole solution (6.11) will be of 
interest. For example, one can easily embed Kerr-Ncwman black hole in this rotating 
monopole space to study a different physical nature of the black holes. 

6.2.2 Kerr-Newman solution 

We can choose the Wang-Wu functions qn{u) such that 

!m, when n — 

-e^/2, when n = -1 (6.12) 
0, when 0,-1 

where m and e are constants. Then we obtain the function from (6.8) 

M(m, r) = m — e^/2r, /x* = a; = 0, 



e2 



and the Weyl scalar is 



which are the same as given (5.1). 



6.2.3 Rotating Vaidya-Bonnor solution 



The rotating Vaidya-Bonnor solution presented above, can also be obtained from 
these rotating Wang-Wu solutions if we choose the functions as 



(ln{u) 



f{u), when n — 

-h{uf/2, when n = -1 (6.14) 
0, when n 7^ 0, —1. 
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Then the corresponding quantities are 



M{u,r) = f{u)-h\u)/2r 



P* = - ^^^2^2 {^^(^/N>" ~ + Q^sin^^ {rf{u),u - hhu),u}, 

This solution includes the rotating Vaidya solution {h{u) — 0) obtained above in (5.9). 
These subsections (6.2.2) and (6.2.3), which are the repetition of the section (5.1) and 
(5.3) above, show that the general rotating solutions with Wang-Wu functions (6.8) 
can employ to derive the rotating Kerr-Newman as well as rotating Vaidya-Bonnor 
solutions in a much simpler way. 

6.2.4 Kerr-Newman- Vaidya solution 

Wang and Wu [32] could combine the three non-rotating solutions, namely monopole, 
de-Sitter and charged Vaidya solution to obtain a new solution representing the non- 
rotating monopole-de Sitter- Vaidya charged solutions. In the same way, we shall 
combine the Kerr-Newman solution with the rotating Vaidya solution obtained above 
in (5.9), if the Wang-Wu functions qn{u) are chosen such that 



Qniu) 



m + f{u), when n — Q 

-eV2, when n = -1 (6.17) 
0, when n ^ 0, —1, 



where m and e are constants and f{u) is the mass function of rotating Vaidya solution 
(5.9). Thus, we obtain the function from (6.8) 

M{u, r) = m + f{u) - e^/2r 

and using this in (6.9) we obtain other quantities 



— i a sin 6* ^ , , 

OJ^—, (6-19 

y/2KRR^-'^ ^ ' 
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The Weyl scalars for this rotating solution are 

This represents a rotating non-stationary Kerr- Newman- Vaidya solution with the line 
element 

ds'^ = [1 - R~'^{2r{m + /(«)) - e^}] du"^ + 2du dr 

+2aR-^{2r{m + /(«)) - e^} sin^^ rfw rf0 - 2a sin^O dr dcf) 

-R^de^ - {(r^ + a^f - A*a2 sin^^} i^-^sin^^ d^^, (6.21) 

where A* = — 2r{m -|- /(«)} -|- -|- e^. Here m and e are the mass and the 
charge of Kerr- Newman sohition, a is the rotation per unit mass and f{u) represents 
the mass function of rotating Vaidya null radiating fluid. The solution (6.21) will 
describe a black hole if m + f(u) > + with external event horizon at r_|_ = 
{m + f{u)} + + /(m)}^ — — e^], an internal Cauchy horizon at r_ = {ni + 

f{u)} — \/[{m + /(m)}^ — — e^] and the non-stationary limit surface r = re{u, 6) = 
{m + f{u)} + \/[{m + f{u)y — a^cos^O — e^]. When we set f{u) = 0, the metric 
(6.21) recovers the standard Kerr-Newman black hole, and if m = 0, then it is the 
rotating charged Vaidya null radiating black hole (5.9). 

In this rotating solution, the Vaidya null fluid is interacting with the non-null 
electromagnetic field whose Maxwell scalar 0i can be obtained from (6.18). Thus, we 
can write the total energy momentum tensor (EMT) for the rotating solution (6.21) 
as follows: 

T„, = T„^)+T„f , (6.22) 

where the EMTs for the rotating null fluid as well as that of the electromagnetic field 
are given respectively 

T^^^ =fi*iah + 2u i^a m^) + 2 a; m,) , (6.23) 



T^^^ = 2p* {£(„n,) + m(,m,)}. (6.24) 

The appearance of non-vanishing uo shows the null fiuid is rotating as the expression 
(6.19) of uo involves the rotating constant a coupling with df{u)/du - both are non- 
zero quantities for a rotating Vaidya null radiating universe (5.9). 
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This Kerr- Newman- Vaidya metric (6.21) can be written in Kerr-Schild ansatz on 
the Kerr-Newman background as 

gT" = + 2g(«, r, ^)44 (6.25) 

where 

Q{u,r,e)^-rf{u)R-\ (6.26) 

and the vector is a geodesic, shear free, expanding as well as rotating null vector of 
both as well as g^^^ and given in (2.8) and is the Kerr-Newman metric (5.4) 
with m = e = constant. This null vector ^a is one of the double repeated principal 
null vectors of the Weyl tensor of g™ . This completes the proof of theorem 1 stated 
above. 

It appears from (6.25) that the Kerr-Newman geometry can be thought of joining 
smoothly to the rotating Vaidya geometry at its null radiative boundary, as shown by 
Glass and Krisch [23] in the case of Schwarzschild geometry joining to the non-rotating 
Vaidya space-time. The Kerr-Schild form (6.25) will recover that of Xanthopoulos 
[33] g'ah = 9ab + when Q{u, r,d) — > 1/2 and that of Glass and Krisch [23] g'^^ — 
g^^^—{2f{u)/r}£aib when e = a = for non-rotating Schwarzschild background space. 
Thus, it can be regarded that the Kerr-Schild form presented in (6.25) above will be 
the extension of those of Xanthopoulos as well as Glass and Krisch. When we set 
a = 0, this Kerr-Newman- Vaidya solution (6.21) will recover to non-rotating Reissner- 
Nordstrom- Vaidya solution with the Kerr-Schild form g^^^ — g^^ — {2f{u)/r}£aib: 
which is still a generahzation of Xanthopoulos and Glass and Krisch in the charged 
Reissner-Nordstrom solution. It is worth to mention that the new solution (6.21) 
cannot be considered as a bimetric theory as g™^ ^ ^{9™ + 9ab)- 

To interpret the Kerr-Newman- Vaidya solution as a black hole during the early 
inflationary phase of rotating Vaidya null radiating universe i.e., the Kerr-Newman 
black hole embedded in rotating Vaidya null radiating background space, we can also 
write the Kerr-Schild form (6.25) as 

9ar = gl + 2Q(r, e)£ah (6.27) 

where 

Q(r, 9) = -{rm - e'^/2)R-^, (6.28) 

Here, the constants m and e are the mass and the charge of Kerr-Newman black hole, 
g^jj is the rotating Vaidya null radiating black hole (5.9) and ia is the geodesic null 
vector given in (2.8) for both g™^ and g^^,. When we set f{u) ~ a — 0, g^fj will 
recover the flat metric, then g^^^ becomes the original Kerr-Schild form written in 
spherical symmetric flat background. 
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These two Kerr-Schild forms (6.25) and (6.27) certainly confirm that the metric 
g™"^ is a solution of Einstein's field equations since the background rotating metrics 
g™ and are solutions of Einstein's equations. They both have different stress- 
energy tensors T^^f-* and T^^^-* given in (6.24) and (6.23) respectively. Looking at the 
Kerr-Schild form (6.27), the Kerr-Newman-Vaidya black hole can be treated as a 
generalization of Kcrr-Ncwman black hole by incorporating Visser's suggestion [34] 
that Kerr-Newman black hole embedded in an axisymmetnc cloud of matter would 
he of interest. Hawking et al [35] have also mentioned the possibility to embed the 
rotating black hole solutions with a theory for which they know the corresponding 
conformal field theory. 

6.2.5 Kerr-Newman- Vaidya-Bonnor solution 

Similarly, one can combine the rotating Vaidya-Bonnor solution obtained above in 
(5.13) with the Kerr-Newman solution (5.4) to generate another rotating solution 
with the mass function 

M{u, r)=m + f{u) - (e^ + h^{u))/2r, (6.29) 

representing a Kerr- Newman- Vaidya-Bonnor solution: 

ds'^ = [l-R-^{2r{m + f{u))-e'^-h^{u)}]du^ + 2dudr 

+2aR-^{2r{m + f{u)) - - h'^iu)} siii^e du dcf) - 2a s\v?e dr dcf) 
-E?de'' - {(r' + a^f - A*a' sin^^} i^-^sin'^ d^^, (6.30) 

where A* = — 2r{m + f{u)} + + + h'^iu). This rotating solution can also be 
written in Kerr-Schild form (6.25) with the function: 

Q{u,r,e)^-{rf{u)-h''{u)/2}R-\ 

When the charge e of the Kcrr-Ncwman solution vanishes, this rotatmg solution (6.30) 
will reduce to a rotating Kerr- Vaidya-Bonnor solution with the mass function: 

M{u,r)^m + f{u)-h?{u)/2r. (6.31) 

It suggests that by choosing the Wang-Wu functions qn{u) properly one can gen- 
erate as many rotating solutions as required. However, the generation of these types 
of rotating solutions will be restricted that the energy-momentum tensor of the fluids 
must be of the form given in (4.1). 
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6.2.6 Kerr-Newman-de Sitter metrics 

Here we shall present the rotating de Sitter as well as Kerr-Newman-de Sitter metrics 
in NP formalism. 

A. Rotating de Sitter solution: 

First we will derive a rotating de Sitter solution of Einstein's equations. For this 
we choose the Wang-Wu functions as 

. . / A76, when n = 3 . . 

= I 0, when n ^ 3 ^^'^^^ 

to obtain the mass function 

A* 

M{u,r) = —-, (6.33) 
o 

The line element for the rotating de Sitter metric is 

A* r'^ 

ds^ = |l -\dv? + 2dudr 

A* r^ 

+2a^- sin^6l du dcf) - 2a sin^O dr d(f) 

-RV - {(r^ + a'^f - A^a^sin^e} R-\m^ed(t)\ (6.34) 

where R^ = r'^ + a^cos^^, A* = — A* r"'/3 + . This corresponds to the rotating de 
Sitter solution for A* > 0, and to the anti-de Sitter solution for A* < 0. In general 
A* denotes the cosmological constant of the de Sitter space. Then the changed NP 
quantities are 



;i-^AV)ri? + A*|, 



0ii = -^^^AVVcos2^, (6.35) 

V'2 = -1 AV^a^cos^g, (6.36) 
A*r'^ 

A = (6.37) 

This means that in rotating de Sitter cosmological universe, the A* is couphng with 
the rotational parameter a. From these NP quantities we can clearly observe that the 
rotating de Sitter cosmological metric is a Petrov type D gravitational field whose 
one of the repeated principal null vectors, is geodesic, shear free, expanding as well 
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as non-zero twist. The rotating cosmological space possesses an energy- momentum 
tensor 



Tab = '2p*i{anb) + 2pm(^amb), 



(6.38) 



where K p* — 2 0ii -|- 6A and Kp = 2 0ii — 6A are related to the density and 
the pressure of the cosmological matter which is, however not a perfect fluid. If 
we set the rotational parameter a = 0, we will recover the non-rotating dc Sitter 
metric [36] , which is a solution of the Einstein's equations for an empty space with 
A = QabR"''^ = A*/6 or constant curvature. However, it is observed that the rotating 
de Sitter metric (6.34) is neither empty nor constant curvature. It certainly describes 
a stationary rotating spherical symmetric solution representing Petrov type D space- 
time. So it is noted that to the best of the present author's knowledge, this rotating 
de Sitter metric has not been seen discussed before. 

B. Kerr-Newman-de Sitter solution: 



By choosing the Wang-Wu function as 
qn{u) = 



m, when n = 

— e^/2, when n = —1 

A*/6, when n = 3 

0, when 0,-1,3 



(6.39) 



we can obtain the function from (6.^ 



,2 ^* ^3 



M(u, r) — m h , 



(6.40) 



where m and e are constants and are the mass and the charge of the Kerr-Newman 
solution. The line element with the function (6.39) is 



ds'^ = {l-R-'^(2mr-e'^ + —-)}du^ + 2dudr 



A* r^- 



+2aR~^ (2mr - H sin^O du dcf) - 2a sui^O dr d(j) 

-RV - {(r^ + aY - A^a^sin^^} R-^ sin^ 6 d(l>^ 



(6.41) 



where ^ + a^cos^O, A* = - 2mr - A* r^/3 + + e^. Then the changed NP 
quantities are 



7 

011 



-=— \(r-m- -A*r^)R - A* 
2RR^ 3 ^ 

l-(e^-AVa^cos^^), 



2R^ R^ 
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A = (6.42) 

We have seen from the above that in each expression of (pn and ip2, the cosmological 
constant A* is couphng with the rotational parameter a. This means that the cosmo- 
logical parameter A* has the effect of its presence in the curvature of the embedded 
Kerr-Newman black hole. The metric (6.41) admits the following energy momentum 
tensor 

Tab = 2p* nb) + 2p m^ojrih) , 
with the density and the pressure of the matter field 



respectively. Without loss of generahty, we can write this Tab with the decomposition 
of p* = + and p = p(^) + that 

Tab = 4p*('^){Aa^f) + HP*^""^ ^(aflb) + 2^^^^ m(„mb)}, (6.43) 

where 



KR^ i?2' 

«*{C) _ ^ „(C) lii_f^2 , ga^cos^^^ 

The advantage of writing Tab in the form (6.43) is that, for Reissner-Nordstrom-de 
Sitter metric (a = 0), the energy- momentum tensor can be written in the form of 
Guth's modification of Tab for inflationary universe [37] as 

Tab^T^^^ +h*gab (6.44) 

where T^^^ is the energy- momentum tensor for non-null electromagnetic field and Qab 
is the Reissner-Nordstrom-de Sitter metric. This indicates that Guth's modification 
of Tab is acceptable only in the case of non-rotating metrics, and its extension in the 
case of the rotating solutions (a ^ 0) will take the form given in (6.43) above where 
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The metric (6.41) describes a rotating stationary solution and is Petrov type D 
with (■02 7^ 0), whose one of the repeated principal null directions is ia- That is, the 
metric can be written in Kerr-Schild form on the de Sitter background as 



KNdS 
'ab 



g^^ + 2Q{r,9)eah 



(6.45) 



where Q{r,9) — —{rm — e^/2)i?~^, and the vector ia is a geodesic, shear free, ex- 
panding as well as rotating null vector of both g^b^'^^ as well as g^f and given in (2.8) 
and g™ is the Kcrr-Newman metric (5.4) with m = e = constant. This null vector 
ia is one of the double repeated principal null vectors of the Weyl tensor of g™"^^ 
and gab- This completes the proof of theorem 2 stated above. We can also write the 
Kerr-Schild form (6.45) on the Kerr-Newman background as 



where Q{r,9) — — (A*r^/6) i?"^, and A* is the cosmological constant and g™ is the 
Kerr-Newman metric (5.4) with m — e — constant. 

It is quite natural to recover rotating de Sitter (m = 0, a 7^ 0, e = 0) and Reissner- 
Nordstrom-de Sitter (m 7^ 0, a = o, e 7^ 0) metrics from the Kerr-Newman-de Sitter 
solution (m 7^ 0, a 7^ 0, e 7^ 0). It is also worth mentioning that from the Kerr- 
Newman-dc Sitter metric, we can recover a rotating charged de Sitter cosmological 
universe (m = 0, a 7^ 0, e 7^ 0). It is noted that one may find the difference between 
this Kerr-Newman-de Sitter metric (6.41) and that of Mallett [8] used by Koberlin 
[38]. Mallett 's derivation of Kerr-Newman-de Sitter metric is based on the direct 
application of Newman- Janis algorithm to the Reissner-Nordstrom-de Sitter 'seed' 
solution. It is also found that the Kerr-Newman-de Sitter solution (6.41) is different, 
in the terms involving the cosmological constant A*, from the one derived by Carter 
[39] and used by Gibbon and Hawking [40], Khanal [41], Hawking, et al [35], and 
others. 

6.2.7 Kerr-Newman-Vaidya-de Sitter solution 

Here we shall combine the Kerr-Newman-dc Sitter solution (6.41) with the rotating 
Vaidya solution given in (5.9). For this purpose we choose the Wang-Wu functions 
qn{u) as follows 



jKNdS 
ab 



= ^„7 + 2Q(r,e)44, 



(6.46) 



■m + f{u), when n = 



qn{u) 



— e^/2, when n = — 1 
A*/6, when n = 3 



(6.47) 



0, when n 7^ 0, —1, 3, 
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where m and e are constants and f{u) is related with the mass of rotating Vaidya 
solution (5.9). Thus, using (6.47) in (6.8) we have the mass function as 



A* 

M{u, r) = m + f{u) - 2^ + —q- 



and other quantities are 



1 



KR^ i?2 



:e^ + AV), 



1 



P = 



KR^ i?2 



{e2_AV(r2 + 2aW)}, 



(6.48) 



UJ — 



KR^ i?2 
— iasin^ , . 



A = g'^'R, 



^2KRR'^ 

A*r2 



ah 



6i?2' 



(6.49) 



and 011, 012, 022 can be obtained from equations (6.48) with (3.10). The Weyl scalars 
are given below 



V'2 

^3 



1 r A*r2 

- i? {m + f{u)} + — -a^cos^^ 



RRR^ 
—i a sm9 



:{(4r + /?)/(«),„}, 



2^2RRR^ 
= ;^=^^{i?V(^),««-2r/(«),„}. 



2i?i?i?2 ^2 

This represents a Kerr- Newman- Vaidya-de Sitter solution with the line element 

A* 



(6.50) 



1 - R-^[2r{m + f{u)) + — e^}] du^ + 2dM(ir 



3 

A* 



where A* 



+2ai?-2{2r(m + f{u)) + — ^ e^} sin^O dudcf) - 2asm'^9 dr d(j) 

-R^de"" - {(r" + a')' - A*a'sin2^} i^-'sin'^d0^ (6.51) 

r2 — 2r{m + f{u)} — A*r^/3 -|- a2 -|- e2. Here m and e are the mass 
and the charge of Kerr-Newman solution, a is the non-zero rotation parameter per 
unit mass and f{u) represents the mass function of rotating Vaidya null radiating 
fluid. When we set f{u) = 0, the metric (6.51) recovers the Kerr-Newman-de Sitter 
black hole (6.41), and if m = 0, then it is the rotating charged Vaidya null radiating 



30 



black hole (5.13). When one sets A* = 0, this metric will recover the Kerr-Newman- 
Vaidya metric ((6.21). In this rotating solution, the Vaidya null fluid is interacting 
with the non-null electromagnetic field whose Maxwell scalar 0i can be obtained from 
(6.48). Thus, we can write the total energy momentum tensor (EMT) for the rotating 
solution (6.51) as follows: 

= + + (6.52) 

where the EMTs for the rotating null fluid, the electromagnetic field and cosmological 
matter field are given respectively 

^i? = yw* + 2 a; mb) + 2 ^(a mb) , 

Tif = 4p*(^){£(,n6) + m(,m,)}, 

^i? = 2{p*(c)£(„n6) + 2p(c)m(„m6)}, (6.53) 
where /i* and cu are given in (6.49) and 



The appearance of non-vanishing uj shows that the null fluid is rotating as the ex- 
pression of ijj (6.49) involves the rotating parameter a coupling with df{u)/du - both 
are non-zero quantities for a rotating Vaidya null radiating universe. 

This Kerr-Newman-Vaidya-dc Sitter metric (6.51) can be written in Kerr-Schild 
form on the de Sitter background as 

9ar'' = 9t! + r, e)eM (6.55) 

where 

Q{u, r, 9) = -[r{m + /(«)} - ey2]R-\ (6.56) 

and the vector is a geodesic, shear free, expanding as well as rotating null vector of 
both g'^^ as well as gfi^^'^^ and given in (2.8). We can also write this solution (6.51) 
in another Kerr-Schild form on the Kerr-Newman background as 

ga^"' = 9a^ + 2g(«, r, e)ia£b (6.57) 

where Q{u,r,9) — —{rf{u) + A*r^/6}i?~^. These two Kerr-Schild forms (6.55) and 
(6.57) certainly confirm that the metric g™^'^^ is a solution of Einstein's field equa- 
tions since the background rotating metrics g™ and g'^^ are both solutions of Ein- 
stein's field equations. They both have different stress-energy tensors T^f^ and T^^^ 
given in (6.53). 
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6.2.8 Rotating Vaidya-Bonnor-de Sitter solution 



We shall combine the rotating Vaidya-Bonnor solution (5.13) with the rotating de 
Sitter solution obtained above in (6.34), if the Wang-Wu functions qn{u) are chosen 
such that 



qniu) 



f{u), when n = 

—e^{u)/2, when n — —l 

AV6, when n = 3 

0, when n^^O, — 1,3, 



(6.58) 



where f{u) and e{u) are related with the mass and the charge of rotating Vaidya- 
Bonnor solution (5.13). Thus, using this qn{u) in (6.8) we obtain the mass function 



M(ii,r) = /(ii)-^ + -g-, 



(6.59) 



and other quantities are 



l-{e^(.)+AV}, 



P 



KR^ 



(6.60) 



1 



K i?2 i?2 



2 r^{/(ii),„ - Je(it) e{u)^u] + a^sin^6' {/(ii),« - ^e{u) e{u)^u\ J , 



UJ 



—i a sin ^ 
y/2KR^R- 



■{i?/(«),„-2e(«)e(«),„}. 



and 011, 012, 022 can be obtained from equations (6.60) with (3.10). The Weyl scalars 
are given below 



V'2 
V'4 



A * 2 

e\u)-Rf{u) + ^a\os^e 



4r + /(«),„} -4e(M)e(M),„ 



— i a sinO 

sin^^^ \^o< r/ \ / \ / \ "1 

— - — - it {r f iu) u — e.(u)e{u) u\ 

2RRR^R^ [ \ M A« \ J \ ),u^^ 

-2r\r f{u)^u- e{u)e{u)^u] ■ 



(6.61) 



This represents a rotating Vaidya-Bonnor-de Sitter solution with the line element 



ds' 



1 - R-'^{2rf{u) + 



)} du^ + 2dudr 
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A* 

+2aR-^[2rf{u) + — e^{u)} sin^ 6 dudc^ - 2a sin^ 6 dr d(j) 

-R^de^ - {(r' + a^f - A*a2sin2^} R-^inH dct>\ (6.62) 

where A* = — 2rf{u) — A*r^/3 + + e^(w). Here, a is the non-zero rotational 
parameter per unit mass and f{u) represents the mass function of rotating Vaidya 
mill radiating fluid. We can write the total energy momentum tensor (EMT) for the 
rotating solution (6.62) as follows: 

Ta, = T;^^+Tif^+Ti^\ (6.63) 

where the EMTs for the rotating null fluid, the electromagnetic fleld and cosmological 
matter fleld are given respectively 

^ib^ = A^*44 + 2u;^(am6) + 2cZ;£(„m6), (6.64) 
ri? = 4p*(^){£(„nfe)+m(„mft)}, (6.65) 
ri?^ = 2{p*(C)£^^^^^ + 2p(C)m(„m6)}, (6.66) 



where 



KR^ i?2' 

A*^4 A* 2 

If we set a = 0, wc recover the non-rotating Vaidya-Bonnor-de Sitter solution, and 
then the energy-momentum tensor (6.63) will take the form of Guth's modiflcation 
of Tab for inflationary scenario [37] as 

Tab = T^b^ + + ^*9a^ (6.67) 

where is the energy-momentum tensor for non-null electromagnetic fleld and gab 
is the non-rotating Vaidya-Bonnor-de Sitter metric tensor. Prom this, without loss of 
generahty, the EMT (6.63) can be regarded as the extension of Guth's modiflcation 
of energy-momentum tensor in rotating spaces. 

The Vaidya-Bonnor-de Sitter metric can be written in Kerr-Schild form 

^r'' = 9a! + 2Q(^, r, e^aib (6.68) 

where Q{u, r, 9) = —{rf(u) — e'^{u)/2} R^^. Here, g'^^ is the rotating de Sitter metric 
(6.34) and ia is geodesic, shear free, expanding and non-zero twist null vector for 
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both g^^ as well as g^^'^^ and given in (2.8). The above Kerr-Schild form can also be 
written on the rotating Vaidya-Bonnor background given in (5.13). 

^rr' = ^J6^ + 2g(r,^)44 (6.69) 

where Q{r,9) — —{A*r^/6)R~^. These two Kerr-Schild forms (6.68) and (6.69) prove 
the non-stationary version of theorem 2 in the case of rotating Vaidya-Bonnor-de 
Sitter solution. If we set f{u) and e{u) are both constant, this Kerr-Schild form 
(6.68) will be that of Kerr-Newman-de Sitter black hole (6.45). The rotating Vaidya- 
Bonnor-de Sitter metric will describe a non- stationary spherically symmetric solution 
whose Weyl curvature tensor is algebraically special in Petrov classification possessing 
a geodesic, shear free, expanding and non-zero twist null vector given in (2.8). One 
can easily recover a rotating Vaidya-de Sitter metric from this Vaidya-Bonnor-de 
Sitter solution by setting the charge e{u) — 0. If one sets a — 0, e{u) = in (6.62), 
one can also obtain the standard non-rotating Vaidya-de Sitter solution [42] . Ghosh 
and Dadhich [43] have studied the gravitational collapse problem in non-rotating 
Vaidya-de Sitter space by identifying the de Sitter cosmological constant A* with the 
bag constant of the null strange quark fluid. Also if one sets a = in (6.62) one 
can recover the non-rotating Vaidya-Bonnor-de Sitter black hole [44]. It certainly 
indicates that all embedded solutions (6.21), (6.30), (6.41), (6.62) can be derived by 
using Wang-Wu functions (6.8) in the rotating solutions (6.4). 



7 Conclusion 

In this paper, we have calculated NP quantities for a rotating spherically symmet- 
ric metric with three variables. With the help of these NP quantities, we have first 
given examples of rotating solutions like Kerr-Newman, rotating Vaidya and rotating 
Vaidya-Bonnor. Then, with the help of Wang-Wu functions, we come to the unpub- 
lished rotating metrics that we combined them with other rotating solutions in order 
to derive new embedded rotating solutions, and studied the gravitational structure 
of the solutions by observing the nature of the energy-momentum tensors of respec- 
tive spacetime metrics. The embedded rotating solutions have also been expressed 
in terms of Kerr-Schild forms in order to show them as solutions of Einstein's field 
equations. 

We would like to mention that Chandrasekhar [20] has established a relation of 
spin coefficients p, /i, r, tt in the case of an affinely parameterized geodesic vector, 
generating an integral which is constant along the geodesic in a vacuum Petrov type 
D space-time 

P = ^ = L = 1, (7.1) 

p fi 71 T 
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This relation is being derived on the basis of the vacuum Petrov type D space-time 
with ^2 7^ 0, Vo = V'l = V'3 = V'4 = and 0oi = 0O2 = 0io = <t>20 = ^12 = ^21 = 
000 = 022 = 011 = A = 0. However, it has been shown in [22] that the non-vacuum 
Petrov type D spacetimes i.e. Kerr-Newman solution possessing electromagnetic 
field and Kantowski-Sachs metric with dust energy-momentum tensor, satisfy the 
Chandrasekhar's relation (7.1). We will show that the very general metric (2.7) with 
three variables, which is of algebraically special in Petrov classification having non-zero 
"02, V's, "04 and stress-energy tensor (4.1), still satisfies the relation (7.1) as follows 

? = ^ = ^ = ? = (7.2) 
p jjt, T R 

This relation (7.2) shows that all rotating solutions, stationary and non-stationary, 
discussed here satisfy the relation (7.1). Thus, it seems reasonable to refer to the 
relation (7.1) as Chandrasekhar's identity as mentioned by Fermandes and Lun [45]. 
This certainly indicates that the NP spin coefficients (3.3) can be used to extend the 
known vacuum results like the relation (7.1) to the non- vacuum ones. Further, we 
also observe from the energy-momentum tensor (4.1) that the metric (2.7) with three 
variables does not include the perfect fiuid. From the above discussion, it certainly 
indicates that the Newman- Janis algorithm can be used to generate rotating solutions 
as shown in section 6, if the metric function M(it, r) is expressed in terms of Wang- 
Wu functions given in (6.8). However such generated rotating solutions from the 
application of Newman- Janis algorithm have limitations that these rotating solutions 
are not included spacetimes admitting rotating perfect fiuid. To have a spacetime 
admitting a rotating perfect fiuid one has to look for another algorithm rather than 
that of Newman and Janis. 

From the above results presented in this paper, it suggests that rotating spacetime 
geometries must also have rotating matter fields, described by the stress-energy ten- 
sors Tab with non zero rotation parameter a. It means that the matching of a rotating 
non-stationary space-time geometry with the usual non-rotating perfect fiuid will not 
have a reasonably good sense. So one needs to look for a rotating perfect fluid to match 
with a rotating non- stationary spacetime geometry. Some of the rotating solutions 
discussed above include rotating non-stationary solutions, like Kerr-Newman- Vaidya 
black hole, Kerr-Newman- Vaidya-de Sitter black hole, Vaidya-Bonnor-de Sitter black 
hole. They possess rotating non-perfect fluids as shown above by the respective Tab- 
To study the nature of these rotating black holes will certainly be a new area of 
interest in classical General Relativity, since all known black hole theorems, like 'no 
hair theorem', Penrose's theorems are based on stationary black holes, rotating or 
non-rotating. 

It is also found that the technique of Wang and Wu with the functions qn(ti) in 
(6.8) can be used to generate rotating solutions as shown in section 6.2 above. By 
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choosing a suitable Wang-Wu function q{u), we obtain a rotating de Sitter space- 
time model. We can also recover the widely used (i) Schwarzschild-de Sitter solution, 
(ii) Reissner-Nordstrom-de Sitter black hole solution, (iii) Kerr-de Sitter solution, (iv) 
Kcrr-Newman-de Sitter solution for early inflation scenarios from the rotating Vaidya- 
Bonnor-de Sitter solution (6.62). These embedded de Sitter space-times can generate 
by using Wang-Wu functions in rotating solutions given in (6.9). This shows that 
Wang-Wu functions in rotating space-time geometry can be applied to generate Kerr- 
de Sitter and Kerr-Newman-de Sitter solutions in a simple way - not directly applying 
Newman- Janis algorithm to Schwarzschild-de Sitter as well as Reissner-Nordstrom- 
de Sitter 'seed' solutions to derive Kerr-de Sitter as well as Kerr-Newman-de Sitter 
solutions respectively. Thus, the rotating solutions (6.9) with Wang-Wu functions can 
avoid the difficulties suggested by Xu [9]. The definitions of embedded spaces used 
here are in agreement with the one defined by Cai et al. [46]. It is worth mentioning 
that the rotating embedded solution, namely Kerr-Newman-de Sitter solution (6.41) 
is found different from the ones discussed in [8,9,40]. Hence, to the best of the 
author's knowledge, the rotating embedded solutions (6.41), (6.51), (6.62) and other 
reducible solutions from (6.62), and also (6.21), (6.30) have not been seen published 
before. Other non-rotating embedded solutions of Einstein's equations can be found 
in Kramer, et al. [47] (and references there in) and Hodgkinson [48]. 

Looking at these overall rotating solutions derived above one can conclude that 

1. all stationary rotating solutions including (a) Kerr- Newman, (b) rotating monopole, 
(c) rotating de Sitter, (d) Kerr- Newman de Sitter solutions which arc derivable 
from the application of Newman- Janis algorithm, arc Petrov type D and each 
spacetime has the repeated principal null vector which is geodesic, shear 
free, expanding as well as rotating. This completes the proof of theorem 3. 

2. rotating Vaidya (5.9), rotating Vaidya-Bonnor (5.13), Kerr-Newman-Vaidya 
(6.21), Kerr-Newman- Vaidya-Bonnor (6.26), rotating Vaidya-de Sitter (6.62) 
when e{u) = and rotating Vaidya-Bonnor-de Sitter (6.62) are all non- stationary 
spherically symmetric solutions. Their Weyl curvature tensors are algebraically 
special in the Petrov classification with null vector ia given in (2.8). This leads 
the proof of the theorem 4 stated in the introduction. 

The remarkable feature of the analysis of rotating solutions in this paper is that all the 
rotating solutions, stationary Petrov type D and non-stationary algebraically special, 
presented here possess the same null vector ia, which is geodesic, shear free, expanding 
as well as non- zero twist. From the studies of the rotating solutions we found that 
some solutions after making rotation have disturbed their gravitational structure. 
For example, the rotating monopole solution (6.11) possesses the energy-momentum 
tensor with the monopole pressure p, where the monopole constant b couples with 
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the rotating parameter a . Similarly, the rotating de Sitter solution (6.34) becomes 
Petrov type D spacetime metric, where the rotating parameter a is coupled with the 
cosmological constant and so on. We have shown that all the rotating embedded 
solutions presented here can be written in Kerr-Schild forms, showing the extension 
of those of Xanthopoulos [33] and of Glass and Krisch [23] . 
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